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Abstract 



The locality hypothesis in density-functional theory (DFT) states that the functional 
derivative of the Hohenberg-Kohn universal functional can be expressed as a local multi- 
plicative potential function, and this is the basis of DFT and of the successful Kohn-Sham 
model. Nesbet has in several papers [Phys. Rev. A 58, R12 (1998); ibid. A 65, 010502 (2001); 
Adv. Quant. Chem, 43, 1 (2003)] claimed that this hypothesis is in conflict with fundamental 
quantum physics, and as a consequence that the Hohenberg-Kohn theory cannot be generally 
valid. We have in a Comment to the Physical Review [Phys. Rev. A 67, 056501 (2003)] com- 
mented upon these works and recently extended the arguments [Adv. Quant. Chem. 43, 95 
(2003)]. We have shown that there is no such conflict and that the locality hypothesis is inher- 
ently exact. In the present work we have furthermore verifled this numerically by constructing 
a local Kohn-Sham potential for the ls2s'^S state of helium that generates the many-body 
electron density and shown that the corresponding 2s Kohn-Sham orbital eigenvalue agrees 
with the ionization energy to nine digits. Similar result is obtained with the Hartree-Fock 
density. In addition to verifying the locality hypothesis, this confirms the theorem regarding 
the Kohn-Sham eigenvalue of the highest occupied orbital. 

pacs: 02.30Sa, 31.15Ew, 31.15Pf 



1 Introduction 



The locality hypothesis is the corner stone of density-functional theory (DFT) and the basis for 
the widely used Hohcnberg-Kohn-Sham model ^El. The central issue is whether the functional 
derivatives of the density functionals used in DFT are strictly local potential functions, and whether 
a Kohn-Sham model with local potential could lead to exact density for ground states of arbitrary 
electronic systems. 

Nesbet has in a series of papers claimed that the locality hypothesis is in conflict with Fermi- 
Dirac statistics and the exclusion [3-20] principle with the consequence that traditional DFT gen- 
erally does not work for systems with more than two electrons. This result has been criticized by 
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Gal [21], by Holas and March [221 123 1 as well as by us |211- We have also recently discussed the 
matter further in a more comprehensive review |2S]. Nesbet has in unpublished works responded 
to the our comment JS| as well as to that of Holas and March . 

The first atomic DFT model was the Thomas-Fermi model, where the electron cloud is approx- 
imated by a homogeneous electron gas. The energy is then a functional of the electron density, 
p(r), which for noninteracting electrons becomes [2111231211 

Ety[p\ = TtfM + j Arv{r) p{r), (1) 

where 

TtfM^C /dr[p(rr/3 (2) 



is the kinetic-energy functional and v{r) the local, external potential. Minimizing this functional 
with respect to the electron density, leads to the Euler-Lagrange equation 

(3) 

dp{r) 

where /i is the Lagrange multiplier associated with the normalization condition, J dr p{r) =^ N . 

The original TF model is rather crude, and does, for instance, not give rise to any electronic 
shell structure. The question raised by Nesbet in several articles [3 is whether there exists 
'an exact Thomas- Fermi model' for noninteracting electrons, which we interpret to imply that a 
functional Ttf[p] exists such that the energy functional Q would be exact for such systems. Nesbet 
claims that such a model would be in conflict with the exclusion principle - the TF equation JSJ 
with a single Lagrange multiplier could never lead to electronic shell structure, which requires one 
parameter for each shell. In the Appendix we have derived the standard Hartree-Fock equations, 
using a single Lagrange parameter, which invalidates the claim that additional parameters are 
needed to generate shell structure. 

Since the existence of an 'exact TF model' is a direct consequence of the validity of the HK 
theorem, Nesbet concludes that the fundamental HK theory is incomplete. For the same reason 
the basic work of Englisch and Englisch |29l nU] as well as most of the current literature |31[ I32L 
ISniiniESlEniEZl is refuted. 

There exists today an overwhelming amount of evidence for the validity of the locality hy- 
pothesis in DFT, and it might seem superfluous to produce another piece. Since the matter is 
still under debate in the literature, however, we shall here try to present simple arguments how 
the conflict has arisen and how it can be resolved. We shall also present in our opinion strong 
numerical evidence in favor of the hypothesis. 



2 Functional derivatives 

Generally, a density functional is defined as a mapping of a normed space of densities^ (M) 

on the space of real numbers {TZ) [55|. 

F[p] : M n, (4) 

a mapping that has to be unique in the sense that a certain density corresponds to a single real 
number. 

^The term 'normed space' implies that all densities of the space have a definite, finite norm, such as ||p{t')|| = 
J dr |p(r)|, but the densities need not be normalized. 
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2.1 The Gateaux derivative 



The Gateaux or weak differential of a density functional F[p] is defined I^S", p. 46], PHI P-293], |37l 
Eq.60] 

,pr .1 F[po + XSp] - F[po] 

dF[po, dp] = lim , (5) 

A — *U A 

or alternatively 

SF[po, XSp] = F[po + XSp] - F[po] = XdF[po, Sp] + uj{po, XSp), (6) 

where 

lim = 0. (7) 

In principle, this differential depends on the direction Sp, and it need neither be linear nor contin- 
uous in Sp. Generally, it can be expressed 

dF[po,M- J drF'{[po,Sp],r)Sp{r), (8) 

where F'(Jpo, 5p],r) is a function of r - not necessarily continuous - and a functional of p as well 
as of Sp. If the differential is linear and continuous in Sp, then the function must be independent 
of Sp, and the differential can be expressed 

dF[po,Sp] = J drF'{[po],r) Sp{r). (9) 

F'([po],f) is here a local, finite, multiplicative function of r that is single- valued for given po, 
and conventionally referred to as the Gateaux derivative of the functional F at the density po |38[ 
p.46], 1^, Eq.61], 

(^) -^'(N,r). (10) 
2.2 The Prechet derivative 

The Frechet or strong differential dF[po,Sp] can be defined p. 37], ISHl p. 292] 

SF[po, Sp] = F[po + Sp] - F[pn] = dF[pn, Sp] + u{pq, Sp), (11) 

where 

lim ^^^^ = 0. (12) 

Here, the limit has to be uniform in the neighborhood of the density po- The Frechet differential 
can also be expressed in the form lO with the derivative in the form of a local potential function (381 
p.41]. 

The criterium for Frechet differentiability H12|l is stronger than the corresponding Gateaux 
criterium jT)) . If the Gateaux derivative exists in the neighborhood of a density and is uniformly 
continuous in this neighborhood, then the Gateaux derivative at that density is also a Frechet 
derivative [38. p. 47], |i39i p. 295]. This requires that the functional is defined also for unnormalized 
densities. In DFT it is possible to stay within the normalization domain, and then Gateaux 
differentiability is sufficient. 

In the next section we shall summarize the basics of the Kohn-Sham model as a background of 
the locality analysis. (For further details regarding the basic theory, see, e.g., refs. |31ll32lll?11l37| .') 
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3 The Hohenberg-Kohn-Sham model 



3.1 The Hohenberg-Kohn theorem 

The Hohenberg-Kohn (HK) universal functional is in the constrained-search formulation given 
by SHEl ^ 

FhkM = mhi{^\f + W\^) = ThkW + WhM (13) 

Here, T and W represent the kinetic-energy and electron-electron-interaction operators of the 
A'^-electron system, respectively. 



N N 

f = -^iv2.; w^J2urzr^r (14) 

i=l i<j ' * -J' 

The normalized wavefunction ^E" belongs to the Sobolev space Hl{Ti?^) ^2 IS ESI EZ], and the 
corresponding functional is defined for all N-representable densities \IV2\ . 

The HK energy functional of a system with an external, local potential v{r) is 

EM = Fhk[p] + j drp{r)v{r), (15) 

and the energy, Eyo, and the electron density, pair), of the ground-state (or the lowest eigenstate 
of certain symmetry 021) are obtained by minimizing this functional over normalized densities 

E^Q = min E^ = i?i,[/Oo]- (16) 

As in our previous works |24l I25j . we shall extend the definition of the HK functional in 
a straightforward way by allowing the wavefunctions to vary also outside the normalized domain 
of the Sobolev space, which is needed in order to be able to apply the Euler-Lagrange procedure. 
The minimization then leads to the Euler-Lagrange equation 

where /i is the Lagrange parameter for the normalization constraint, J dr p{r) — N . 

If we omit the electron- interaction part, W , of the HK functional then the EL equation 
(|17|l is quite analogous to the Thomas- Fermi equation for noninteracting electrons ■ The kinetic 
part, Thk[p], of the HK functional (|13|l is then the kinetic-energy functional of the 'exact Thomas- 
Fermi model', which thus is a direct consequence of the HK theorem. 

3.2 The Kohn-Sham model 

In the Kohn-Sham (KS) model, the interacting system is replaced by a system of noninteracting 
electrons, moving in the local Kohn-Sham potential, VYis{i')^ 

[-\V2 + VKsir)] Mr) = <^,(r). (18) 

The analogue of the energy functional H15|l is then 

E,[p] = Ts[p] + J drpir) VKsir), (19) 



^This is known as the Levy-Lieb functional and sometimes denoted -Fll[p]- 
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where 

TsM = miii($|f|$), (20) 

and $ is a single Slater-determinantal wavefunction in the same functional space as previously. 
Minimizing the energy functional, leads in analogy with the Euler-Lagrange equation (|17|l to 

S + -Ks(r) = /.. (21) 
The two equations lead to the same solution, if - apart from a constant - 

If this potential were known, solving the KS equations self-consistently would, in principle, yield 
the exact electron density and, if the HK functional 115|) were known, the exact total energy of the 
ground state - or the lowest state of a given symmetry. In addition, it has been shown that the 
eigenvalue of the highest occupied KS orbital yields the exact ionization energy |43[ I44[BS1|36| . 

A crucial point in the KS procedure is that the functionals involved are differentiable and that 
the functional derivative is in the form of a local potential function. This is now a well-established 
fact. In the standard DFT procedure, the variations are restricted to normalized densities, and 
then Gateaux differentiability is sufficient. The Gateaux differentiability of the functional 1)1 3|l was 
rigorously demonstrated two decades ago by Englisch and Englisch [211 EOI , based upon works of 
Lieb |34| . and the derivative was shown to be in the form of a local potential function.'^ These 
results have recently been carefully confirmed by van Leeuwen [37| . By extending the definitions to 
unnormalized densities, as mentioned above, we have demonstrated that the functionals involved 
have a Gateaux derivative according to the definition Q above also in the extended domain. This 
derivative is - at least for practical purposes - also a Frechet derivative [2^123 EE] 



4 Nesbet's locality dilemma and its solution 
4.1 The problem 

In his efforts to demonstrate that standard DFT is in conflict with fundamental physics, Nesbet 
considers the ground state - or lowest state of given symmetry - of a system of noninteracting 
electrons. For a two-electron system this can be represented by a determinantal wavefunction 
^(''i: '"2), composed of spin-orbitals with the space part satisfying the single-electron Schrodinger 
equation 

(i-f f(r))(/.,(r) =£,;0,(r) (23) 

with i being the kinetic-energy operator i = — ^V^ and v{r) a local external potential. For the 
kinetic energy and electron density of the system Nesbet uses the expressions 

T = (0i|£|0i) + (02|i>2> (24) 
p{r) = p^{r)+p,ir) = \Mr)\^ + \Mr)\', (25) 

valid for normalized orbitals. By making small orbital changes the kinetic energy of the ground 
state is to leading order changed by (Eq. 9 in [2]) 

dT = j dr{ei-v{r))5pi{r) + J dr {e2 ~ v{r)) Sp2{r), (26) 

^The proof concerns the so-called Lieb functional, usually denoted -Fl[p]i which in the case of nondegenerate 
ground states we are concerned with here is identical to the functional IlijI . 

''This holds for calculations with finite basis set, where the kinetic energy has a finite upper bound, but it may 
not be rigorously true when the number of dimensions is infinite. (See ref. l'2bl for more detailed discussion). 
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where 

Spi{r) = 5(j)*{r)(j)^{r) + c.c. 

This is neither of the form ^ nor of the form © outside the domain of normahzed orbitals 
(where J dr5pi{r) ^ 0) when the eigenvalues are different. Therefore, dT cannot be a differential 
of a density functional in this region, and no density-functional derivative exists, as we have 
demonstrated before |24|. Instead, Nesbet introduces 'orbital- dependent' derivatives 

- v{r) (27) 



and interprets them as Gateaux derivatives, which is obviously neither in accordance with the 
standard definition (jS)) above, nor with the extended definition JSJ. 

In spite of the fact that T is not a density functional in the extended domain, Nesbet applies 
the chain rule, 

5T ^ JT_ 5p{rl ^ JJ_^ 
5pr{r) 5p{r) Sp,{r) 5p{r)' 

which together with the result (|27|l leads to a conflict with the TF equation JJl, if the orbital 
energies are different, 

ST , ST 
Spi{r) dp2{r) 

The reason for the conflict is that the kinetic-energy expression H24I) . which obviously is an or- 
bital functional, is not a density functional outside the domain of normalized orbitals when the 
eigenvalues are different. Therefore, the chain rule 1)281) cannot be used in that region, and the 
conflict is only apparent, as we have demonstrated before and emphasized also by Holas 

and March [221 ■ shall now indicate how this conflict can be resolved. 



4.2 The solution 

In the treatment of Nesbet the density variations are allowed to go outside the normalization 
domain. At the same time it is emphasized that the kinetic-energy expression should be an orbital 
as well as a density functional, which, of course, has to be the case in the extended domain. 

The kinetic energy and the electron density are not a priori defined outside the normalization 
domain, and we have to choose a proper extension. A natural extension is simply to use the 
standard expressions also outside the normahzation pHI^ . 



JJ dndrs **(r-i,r2)(ti +t2)$(ri,r2) = {(j^M + {^2^2} (0i|0i) (30) 

p(r) = 2 J dr 2 Mr, r2)\^ ^\Mr)\' {c^M + \Mr)m<l>i\<Pi), (31) 

assuming the orbitals to be mutually orthogonal. Since the orbitals are not normalized, it is vital 
to maintain the normalization integrals in the expression to be varied. Instead of the expression 
(I26|l . the differential of the kinetic energy due to orbital changes now becomes [2] 

dT= [dr (^i±^ - v{r)) Sp{r), (32) 



which is in accordance with the expression lO also for variations outside the normalization. This 
verifies that the form is a unique functional of the total density (|31)) in the neighborhood of 
the ground-state density, and it yields the derivative 

ST ei+e2 , . , . 

— o v{r). (33) 



5p{r 



This is in contrast to the result of Nesbet H29() orbital independent. The derivative is the same in 
all directions - also outside the normalization - and is, as discussed above, of Frechet type. 
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We can now identify the constant in the derivative with the Lagrange multipUer in the 
'Thomas- Fermi' equation 

M=^— ■ (34) 

Thus, there is only a single Lagrange parameter, when the extension to the unnornialized domain 
is done in a consistent way, and the conflict (|29|) with the Thomas-Fermi equation is removed. 

The crucial point in the arguments is the chain rule (|28|l . The more exact form of this rule is 



ST 



dr' 



5T 5p{r') 



5p{r') 5(1)* {r)' 



(35) 



If we assume that the density is given by then = 5{r-r') (j)i{r) and = 4>i{'r), 

which is equivalent to the simplified chain rule H28|l . But this holds only inside the normalization 
domain, where the density can be written in the form (|25|l . In the treatment of Nesbet, however, 
the variations have to go outside this domain, otherwise the orbital dependence of the differential 
(|26|l would vanish, and there would be no conflict. 

With the expression for the density, valid also outside the normalization domain, the 
correct chain rule 1)35(1 yields 



6T 



ST 
Sp{r) 



(</'2|'^2) + U 



ST 



This is compared with the direct orbital derivation |24| 
ST 



Sp 



Mr)- 



El + £2 



«(r))(02|02)-(02k|</'2) <i^i{r), 



(36) 



(37) 



and identification yields the same expression for the derivative as before . Thus, by using proper 
expressions and using the chain rule in a proper way, the conflict observed by Nesbet between direct 
orbital derivation and the chain rule has disappeared. 



4.3 Comments 

The basic claim of Nesbet is that 'density functional derivatives of the Hohenberg-Kohn universal 
functional cannot be equivalent to local potential functions, as assumed in established literature 
l.'yjjf . for more than two electrons in a compact electronic system' |15j . and it is concluded that 'the 
TF theory is not equivalent to Kohn-Sham theory, even if the exact Hohenberg-Kohn universal 
functional were known and used' ;3 ■ 

Reference is made to the 'well-known failure' of Thomas-Fermi theory to describe the atomic 
shell structure. This failure is unknown to us - unless, of course, it refers to the trivial case of the 
original approximation (0) from the 1920's. According to the HK theorem, an exact kinetic-energy 
functional does exist (ThkM in Ea. (|13() ). but to our knowledge no approximation to the functional 
beyond the original approximation has been constructed and tested. 

A related question addressed by Nesbet is that the Euler-Lagrange equation, like the TF equa- 
tion ((2Jl or the HK equation IjlTfl , with a single Lagrange parameter cannot lead to electronic shell 
structure. In a recent paper 20, it is stated that 'the exclusion principle requires independent 
normalization of the orbital partial densities'. Without such independent normalization Nesbet 
claims that the electronic wavefunction will collapse to the lowest single-electron state. 

In the Appendix we have given an elementary counterargument by deriving the standard 
Hartree-Fock equations for a two-electron system with a single parameter for the normalization of 
the total wave function. This obviously leads to shell structure, and there is no collapse into the 
lowest electronic state 18 . Parameters for normalization of individual shells are not needed, if 
orbital overlap integrals are included in the expression to be minimized. 
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The standard way of deriving the HF equations is, of course, to introduce Lagrange muhiph- 
ers for the orthonormahty condition. Here, however, we want to emphasize that there exist two 
equivalent ways of performing a minimization by means of orbital variations, either by enforcing 
orthonormality by means of Lagrange multipliers, omitting orbital normalization and overlap inte- 
grals, or by maintaining these integrals together with a single Lagrange parameter for the overall 
normalization (and performing the orthonormalization afterwards, if desired). In our DFT work 
we used the latter method. 

The reason for the failure of standard DFT is according to Nesbet that the treatment is re- 
stricted to normahzed densities ^lEl- The Euler-Lagrange procedure requires the functional to 
be defined also in infinitesimal neighborhoods of such densities. This is regarded as 'a crucial lim- 
itation of the Hohenberg-Kohn theory. ' It is assumed that the normalization constraint 'hides' the 
additional parameters needed to generate shell structure. The conclusion is that a correct theory 
can only be developed using functional derivatives of orbital densities with nonlocal potentials. 

We have shown |24[ 05] that restricting the treatment to normalized densities constitutes no 
limitation. Extending the treatment to the unnormalized domain leads to identical results - no 
'hidden' parameters are uncovered. 

That DFT is complete without any additional parameters has also been demonstrated among 
others by Perdew and Levy ^HI, who emphasize that 'the exact Junctionals themselves (and not 
extraneous constraints ) are responsible for the shell structure '. 

Our analysis reconfirms the well-established result that DFT with local potential is inherently 
exact. The same conclusion, using different lines of arguments, is reached by Holas and March in 
their recent analysis 



5 Numerical demonstration of locality 

As a further demonstration of the validity of the locality hypothesis, we have, using the method 
of Baerends and van Leeuwen |47| . constructed numerically the local Kohn-Sham potential for 
the 1525^5 state of neutral hehum - a state often used by Nesbet to demonstrate the break- 
down of this principle. Our starting point is the density generated by an accurate many-body 
wavefunction, obtained by means of our all-order pair procedure The potential is constructed 
so that the Kohn-Sham orbitals generate the many-body density very accurately. The energy 
eigenvalue of the 2s Kohn-Sham orbital is then found to agree with the many-body ionization 
energy within our numerical accuracy of nine digits |49| . This result verifies that a single local 
potential can generate essentially exact results, also for electrons with different energy eigenvalues, 
and that there is no conflict between the locality hypothesis and the exclusion principle, as claimed 
by Nesbet. Related results (with less numerical accuracy) have been obtained by many groups 
in the past mEOlEIlESESlElEHlESIEII- In addition to verifying the locality hypothesis, 
our result represents another numerical verification of the theorem regarding the highest occupied 
Kohn-Sham eigenvalue 03 El ESI EH EHI • 

We have obtained similar result with the Hartree-Fock density, which also can be reproduced 
by a local Kohn-Sham potential. The eigenvalue of the KS 2s orbital then agrees with the HF 
value with high accuracy. This demonstrates that it is possible to represent exchange - in the DFT 
sense - by a local potential. 

The latter result is in sharp contrast to the conclusions drawn, for instance, by Nesbet and 
Colle who conclude that the Kohn-Sham model with a local potential can never be 

exact and, in particular, that a local potential can never be better than the optimized-potential 
method of Talman and Shadwick (SHI- This potential in turn yields higher ground-state energy 
than the Hartree-Fock model, as measured by the expectation value of the Hamiltonian, {H), and 
this is supposed to confirm the failure of the locality hypothesis for exchange. A comparison of 
this kind is, of course, of no value, since the energy in DFT is given by the HK energy functional 
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- not by {H). Needless to say, though, that a local potential can never in all respects reproduce 
the nonlocal exchange. 



Our results provide counterexamples to Nesbet's conclusions drawn in several papers |14l 
that a correct theory can only be developed using or&zto/-functional theory with nonlocal potentials. 
We have demonstrated that density- functional theory with local potential is exact - but, of course, 
not identical to orbital-functional theory with nonlocal potential. 



6 Summary 

In summary, we have reconfirmed that the locality hypothesis is an exact statement. The conflict 
with the exclusion principle, claimed by Nesbet, is only apparent and is eliminated by using ex- 
pressions that are density functionals also outside the orbital normalization. Our conclusion is 
supported by accurate numerical results. 
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Appendix: Fermi-Dirac Statistics with a single Lagrange Mul- 
tiplier 

In this appendix we want to demonstrate in an elementary way that it is possible to satisfy 
the conditions for Fermi-Dirac statistics in a variational procedure with only a single Lagrange 
multiplier. We shall do that by minimizing the expectation value of the Hamiltonian {'i>\H\'^) 
under the condition that the total wavefunction "if is normalized, (^'l^') = 1, and show that this 
leads to the standard Hartree-Fock equations. No additional parameters are needed. Fermi-Dirac 
statistics is enforced by demanding that the two-electron wavefunction is antisymmetric. 

Generally, the solution of the problem is found by demanding that 



is stationary with E being the Lagrange multiplier for the normalization of ^E*. This leads to the 
Schrodinger equation 



We consider now a two-electron system in a state where the electrons arc in different orbitals, 
such as the lowest triplet state, ls2s^S. The antisymmetric wavefunction can be expressed in 
terms of spin-orbitals (j)a and (/){,, 



(38) 



(39) 




(40) 



The hamiltonian H is here 



H 



ho{l) + ho{2) + — 



(41a) 



9 



with 

ho = -lw'--. (41b) 
2 r 

Using Eq. (|40|l . the functional H38() now becomes an orbital functional 

i[*HF] =i[</'a,06] = {a\ho\a){b\b) + {a\a){b\ho\b) - {a\ho\b){b\a) - {a\b){b\h„\a) 

+ (ab\ — \ab) - (ab\ — \ba) - E((a\a)(b\b) - (a\b)(b\a))^ (42) 

which should be stationary under independent small variations of the orbitals. So, the conditions 
SL/5(p'^ = and SL/5(f>l = lead to the orbital equations 

ho\a){b\b) + \a){b\ho\b)-ho\b){b\a)-\b){b\ho\a) 

+ {b\±\b)\a)-{b\ — \a)\b)-E(\a){b\b)-\b){b\a))=0 (43a) 

ho\b){a\a) + \b) {a\ho\a) - ho\a) {a\b) ~ \a) {a\ho\b) 

,-|a)|6)-(a|— , 
ri2 ri2 



{a\ — \a)\b) - {a\ — \b)\a) - E(\b){a\a) - \a){a\b)) = 0, (43b) 



while SL/S(j)a — and dL/S(j)b — lead to equivalent complex conjugate equations. 

Although the normalized two-electron wave function is unique (up to a phase factor of modulus 
one), the orbitals are not. They need not be individually normalized and/or orthogonal. If we have 
found orbitals that satisfy the equations above, however, we can easily transform them to become 
orthonormal without changing the two-electron wavefunction. Denoting the transformed orbitals 
again as \a) and \b) for brevity, we insert their properties 

{a\a) = {b\b) = 1 (44a) 
{a\b) = {b\a) = (44b) 

into the equations H43a|l and Ij43bp to obtain 

ho\a) + \a){b\ho\b)-\b){b\ho\a) 

+ — |6)|a)-(6| — |a)|6)-^|a)=0 (45a) 

ho\b) + \b){a\ho\a)~\a){a\ho\b) 

+ (a|^|a)|5)- (a| — |6)|a) =0 (45b) 
ri2 ri2 

which are equivalent to the familiar HF equations 

ho\a) + {b\ — \b)\a) ~ {b\ — \a)\b) = ea\a) + eba\b) (46a) 
ri2 ri2 

ho\b) + {a\ — \a)\b) - {a\ — \b)\a) = eb\b) + eab\a) (46b) 
ri2 ri2 

if we denote some combinations of constants as 

Sa = E-{b\ho\b) (47a) 
Eb = E-{a\hn\a) (47b) 
Sab = {a\ho\b) = e*b^. (47c) 

By projecting (|45a|l onto (a| or (|45b|l onto (&|, the value of the Lagrange multiplier E becomes 

E = {a\ho\a) + {b\ho\b) + {ab\ — \ab) - {ab\ — \ba), (48) 

ri2 
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which shows that E corresponds to the total energy, E — (^'Hrl-ff I^hf)- By inserting this value 
into H47a|l and Ij47b|) we also get the orbital energies Sa, £b 

Ea = (^a\hQ\a) + {ab\ — \ab) - {ab\ — \ba) (49a) 

£b = (blholb) + (abl — lab) - (abl — lba) (49b) 
ri2 ri2 

An alternative form of the total energy (|48|l is therefore 

E = £^ + £^- ({ab\ — \ab) - (ab\ — \ba)). (50) 
V ri2 ri2 ' 

Note that it is not possible to introduce the simplification, i.e. orbitals transformed into or- 
thonormal orbitals, already in Eq. H42|) . before the variation. That would not lead to meaningful 
orbital equations. It is the variation of the overlap integrals that is responsible for making the 
orbital energies £a and e\y different. 
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